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ABSTRACT 


Predictions of propagation loss made by the computer programs Radio Physical Optics (RPO) 
Computer Software Configuration Item (CSCI) and M-Layer are compared. The results of the high 
frequency parabolic equation approximation, as formulated in RPO, agree almost always with 
those derived from the low frequency modal computation as fonnulated in M-Layer. But at low 
altitudes in the neighborhood of the radar horizon, deviations between RPO and M-Layer become 
significant for some cases. RPO appears not to be able to properly account for the effects of a 
high altitude surface-based duct at a short range. Since the discrepancies fall in regions of 
importance to naval operations, a definitive resolution is an urgent task to be undertaken in the 
immediate future. 
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I. INTRODUCTION 


A. THEORETICAL BACKGROUND 

The importance of environmental effects on communication links and radar systems 
has been recognized for well over two decades. The proliferation of computer programs 
such as the Integrated Refraction Effects Prediction System (IREPS) [Ref. 1] and the 
Engineer’s Refractive Effects Prediction System (EREPS) [Ref. 2] within the Navy for 
the prediction of propagation loss of radio waves, together with the general availability 
of computing power brought along by the personal computer (PC) revolution, has greatly 
increased the Navy’s awareness of such effects. 

EREPS essentially is a PC version of IREPS, with provisions for greater flexibility 
in setting input parameters. It employs a combination of different methodologies to 
deduce the propagation factor: ray-optics within line-of-sight; curve-fitting, together with 
frequency scaling based on the assumption that a single mode contributes to the complete 
field strength [Ref. 3J, output from M-Layer [Ref. 4, 5] in the over-the-horizon region, 
and linear interpolation of the results from ray-optics and M-Layer in the region in 
between. Through such curve-fitting techniques, EREPS (and IREPS) gains speed at the 
expense of accuracy. 

The M-Layer program, though time-consuming to run, provides results beyond the 
crude approximation of EREPS, especially in the penumbra region where hundreds of 
modes are required to give an accurate reading. As can be seen in Chapter HI, M-Layer 
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works well at low altitudes, even in the illuminated region. The theory of M-Layer is 
formulated in terms of the excitation and propagation of electromagnetic waves in the 
earth-atmospheric waveguide. The earth, in fact, is treated as a flat, lossy, homogeneous 
half-space. Its curvature is compensated for by postulating the existence of a linearly 
increasing component of the index of refraction in the atmosphere, in addition to the 
slight natural variation of the index of refraction of real air. The rate of increase of the 
linear component is set to the inverse of the effective earth radius. This new index of 
refraction is called the modified index of refraction, and this method of simplifying the 
problem constitutes the earth-flattening approximation [Ref. 6]. Combined with the 
piecewise-linear approximation to the modified index of refraction, the eigenfunctions of 
the flattened earth-atmospheric waveguide are Airy functions of different variables over 
different regions of the atmosphere of linearly varying modified refractive indices. The 
waveguide modes are found numerically, which is the process taking up most of the 
computation time. 

Recently, the RDT&E Division of the Naval Command, Control and Ocean 
Surveillance Center (NRaD), which previously produced IRE PS, EREPS and M-Layer 
when it was known as Naval Ocean Systems Center (NOSQ, came out with a Radio 
Physical Optic (RPO) program [Ref. 7}. It is structured in the same spirit as EREPS: 
different approximations are used over different regions of space. The pan in EREPS 
which relies on curve-fitting of M-Layer predictions is replaced with direct computation 
based on the parabolic equation approximation (PE) [Ref. 8, 9]. This makes RPO far 
superior to EREPS. Computational expediency precludes extensive application of PE. 
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Instead, PE computations are launched at a range of 2.5 km, unless the incident ray to the 
earth gets near the diffraction limit first. A split-step, outward advancing scheme for 
solving the parabolic equation is utilized. At each step, a fast Fourier transform (FFT) is 
performed. The FFT size is set between 7 to 10 powers of 2, while the separation between 
FFT data points is limited to within 2 to 49 wavelengths. Hence, PE computations will 
never be done above a height of 50,176 wavelengths. At 10 GHz, this represents a height 
of no more than 1.50528 km. To avoid aliasing, the FFT coefficients computed from 3/4 
of the height to the top of the region where the FFT is carried out are filtered with a 
factor which decreases from 1 to 3/4 following a sine-squared variation. Furthermore, the 
initial field strength at the range where PE computations start is weighted by a Gaussian 
factor which decays to -70 dB over the same altitudes. Thus, only the lowest 3/4 of the 
FFT results are actually used as PE solutions. 

To fill the remainder of die space, RPO relies on ray-optics until the incident ray 
to the earth approaches the diffraction limit Beyond this limit PE results are utilized 
where they art* available. This is called the PE region. Above the PE region, if the space 
is not covered by ray-optics, an extended optic (XO) region is defined. The PE solutions 
at the greatest height are treated as rays emanating from below, their interferences with 
the direct rays, if present are taken as the fields in this region. 

Even within the ray-optics region, RPO subdivides this region into flat-earth (FE) 
and ray-optics (RO) regions. The earth is considered as flat and all refraction effects are 
ignored in the FE region, which is limited to 2.5 km in range unless the antenna is 
elevated beyond 5°. 
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B. HIGH-FREQUENCY VERSUS LOW-FREQUENCY APPROXIMATION 


Beyond regions where ray-optics is applicable, RPO relies on PE. The parabolic 
equation approximation to the Maxwell wave equations is developed under the optical 
assumption that the operating frequency is so high that a main direction of wave 
propagation in terms of a plane wave can be specified. Hence, only the variation of the 
"envelope" which modulates the magnitude and phase of the plane wave has to be 
considered, instead of the fast variation of the complete wave to the order of the 
wavelength. On the other hand, M-Layer has its theoretical basis in an eigenfunction 
expansion. This is an approach most suitable for low frequency applications because more 
and more modes will be needed as the size of the geometrical structure becomes large 
compared to the wavelength. Even though M-Layer removes the earth radius as the 
dominating length scale through the adoption of a flattened earth, the thickness of the 
overall atmospheric layers specified is still far greater than the wavelength in all 
applications. The comparison of the predictions of PE against M-Layer is one between 
two theories which approach a problem from two extremes. The fact that the results of 
the high frequency PE computations as formulated in RPO agree almost always with those 
derived from the low frequency modal computation as formulated in M-Layer represents 
a significant engineering achievement by NRaD. However, there are disagreements at low 
altitudes in the neighborhood of the radar horizon under some ducting conditions. These 
are regions of importance to naval operations. A definitive resolution of these 
discrepancies is an urgent task to be undertaken in the immediate future. 
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C. SCOPE OF THE WORK 


In this thesis, the radiation of a vertically polarized transmitting antenna located at 
a height of 15 meters is investigated. Three modified refractive index profiles, with the 
first two given in Ref. 7, are utilized: a 300 m surface-based duct, a 14 m evaporation 
duct, and the combination of these two profiles with one sitting on top of the other. Three 
frequencies at 3 GHz, 6 GHZ and 12 GHz are chosen which give an approximate 
coverage of the radar bands. Since only low altitude, near the horizon propagation is of 
interest, propagation loss for altitudes up to 100 meters within the range of 15 to 110 km 
are evaluated by both RPO and M-Layer for comparison. 




a PROGRAM SETUP 


A. TEST CONSIDERATIONS 

Of interest to this thesis research are the propagation loss predictions made by the 
parabolic equation approximation method and by waveguide mode computations. The 
NPS version of M-Layer is a waveguide mode computation program which has been 
greatly enhanced in its accuracy and efficiency [Ref. 5] from its original incarnation 
assembled by NOSC [Ref. 4]. As explained in Chapter I, the RPO program, completed 
by NRaD and available to NPS only recently, includes several methods of computation 
applied to different regions of space. To ensure that only predictions made under the 
parabolic equation approximation are included in the comparison, the RPO FORTRAN 
source code [Ref. 7] is used. This source code consists of a set of subroutines. The output 
for each region comes from a distinct subroutine. By instructing every RPO subroutine 
which computes the propagation loss to associate a code number to its output, the 
predictions made with the parabolic equation approximation are easily identified. 

The RPO FORTRAN source code comes without a main program. Nor does it 
contain input routines and files. Thus, the first task of this research was to write a main 
program for RPO which will call the relevant subroutines to compute the desired output 
Next, data files had to be created to provide the program with refractive index profiles 
and other necessary parameters. An input subroutine had to be written to feed these data 
to the main program and the RPO subroutines. Furthermore, to ensure that both RPO and 






M-Layer use the same environmental data and compatible input parameters, a new M- 
Layer input subroutine was created which replaces the subroutine of the same name in 
the NPS version of M-Layer. Finally, the output data are presented graphically using a 
plotting routine written in MATLAB, version 3. Specific considerations put into these 
subroutines are discussed in the following sections, with the RPO main program 
"RPOmain," the RPO input routine "RPOstdin," the substituting M-Layer input routine 
"MLstdin," listed respectively in Sections 1 through 4 of Appendix A. 

B. RPO MAIN PROGRAM 

In RPO, the parabolic equations are solved with the split-step algorithm. As a 
prerequisite, RPO has to estimate the electric field at an initial range. It then proceeds, 
at a pre-determined step size, toward the range where a propagation loss prediction is 
called for, computing the electric field at every step along the way. Hence the subroutines 
of RPO can be separated into two types: those providing other subroutines with initial 
values and defining their boundaries of application, and those handling the PE 
computation at each step. After the tedious but necessary bookkeeping of declaring global 
variables and setting up common blocks to pass along values, the main program first 
reads in data by calling the input subroutine "RPOstdin” and initializes the subroutines 
by calling "RPOini" before entering into a loop to step the computation through the 
desired range and write the propagation loss predictions within the selected region to an 
output file. A flow chart for the main program, RPOmain, is shown in Figure 1. A listing 
of the program is included in Appendix A. 
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C. INPUT FILES AND SUBROUTINES 


RPO and M-Laycr have different capabilities and require different types of input 
For example, RPO needs a presumed radiation pattern to set up its initial electric field to 
begin computation; it cannot compute electromagnetic (EM) radiation directly from a 
specified current distribution. On the other hand, M-Layer computes the radiation from 
a dipole source; it thus can not simulate antennas arbitrarily. RPO is able to handle 
multiple refiractivity profiles at different ranges due to further approximations made in the 
theory beyond the reduction of the wave equation to a parabolic equation while M-Layer 
is restricted to deal with a single protile. Nevertheless, they use many parameters in 
common. Since the purpose of this work is to compare the predictions of these two 
programs, it is desirable to have them read in exactly the same tiles whenever possible. 
This is especially true for the refiractivity protiles. 

For this thesis, three ASCII files are set up; one contains the parameters common 
to both RPO and M-Layer, including the modified rcffactivity profile; another contains 
RPO specific parameters such as the antenna pattern; the third contains parameters used 
only by M-Layer such as the parameter aloss, which specifies the greatest range 
attenuation rate of the modes to be searched, in dB per kilometer. The parameters which 
are common to both RPO and M-Layer are listed in Table HI. Those specific to RPO 
are listed in Table IL2. The M-Layer specific parameters are listed in Table n.3. Two 
batch files are written to combine automatically the two proper files for input into RPO 


and M-Layer respectively. 





Table 0.1. Input parameters common to RPO and M-Layer. 


RPO 

M-Layer 

Description 

nlevls 

nzlayr 

number of profile levels (a)(b) 

wind 

wind 

surface wind speed 

fmhz 

fqmzin 

operating frequency 

ipolar 

mpol 

antenna polarization^ 5 ^ 

ztran 

ztinit 

transmitter height 

mlxout 

nx 

number of receiver ranges 

xinit 

xinit 

initial range of the receiver 

delx 

delx 

receiver range increment 

nzout 

nzr 

number of receiver heights 

zrinit 

zrinit 

initial height of the receiver 

delzr 

delzr 

receiver height increment 

zprof 

zi 

heights at which profile data 

are specified (an array) 

capm 

zim 

profile data (modified index of 

refraction; an array) 


(a) RPO: from 1 to nlevls; M-Layer from 0 to nzlayr. 

(b) follows M-Layer definition, adjust RPO during input 

(c) RPO: 1 horizontal; 2=vertical; M-Layer O=horizontal; l=vertical. 
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Table n.2. RPO-specific input parameters. 


RPO 

Description 

Remarks 

nxout 

number of output ranges computed 

calculated by Xmax/delx 

where 

Xmax=Xinit+delx*mlxout 

selx 

minimum range to output data 


xprof 

ranges at which index of refraction 

profiles are specified 

set to 0.0 

nprofs 

number of profiles specified 

set to 1 

ipatm 

antenna pattern 

1: omni-direction 2: sin(x)/x 

3: Gaussian 4: cosecant-squared 

5: height-finder 6: user defined 

use 1 or 3 only 

beamw 

antenna elevation beam width 


elang 

antenna elevation angle 

set to 0.0 

nfacs 

number of height-finder data 

set to 0 (not used) 

iscatt 

includes troposcatter 

0: no 1: yes 

set to 0 (not used) 




Table 0.3, M-Layer-specific parameters. 


M-Layer 

Description 

Remarks 

mfile 

0: read input and compute eigenvalues 



1: read eigenvalues as input 


delfq 

frequency increment 

set to 0.0 (not used) 

nfreq 

number of frequencies to be used 

set to 1 

aloss 

maximum range attenuation rate in 

set to 2.0 or 5.0 


dB/km of modes to be found 


debt 

height increment of transmitter 

set to 0.0 (not used) 

nzt 

number of transmitter heights 

set to 1 

refz 

reference height at which refm and 

set to 0.0 


refgab are given 


refm 

modified refractivity at the 

set to 339.0 


reference height 


refgab 

modified gas absorption at the 

set to 0.0 (not used) 


reference height 
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In RPO, special formulas for the dielectric constant and conductivity of sea water 
as a function of frequency are used. Sea surface roughness is also given in terms of a 
function of wind speed. To incorporate this into M-Layer, the input subroutine "MLstdin" 
for M-Layer is modified to carry out these computations using the same formulas before 
providing these parameters to the program, even though surface roughness is not 
considered in this comparison. 

D. VALIDATION 

Before proceeding with the comparison, a few RPO test cases are run using this 
newly written main program and the input subroutine. Specifically, the tests listed in 
Ref. 7 under the names LOBW (low beam width limit), SBDUCT (surface-based duct) 
and EDUC7 (evaporation duct) are carried out. The tabulated results in Ref. 7 are 
reproduced exactly as long as the range increment delx and height increment delz arc the 
same as those specified therein. Otherwise, variations of up to 4% are observed. 



IE. RESULTS 


Propagation of waves through several refractivity profiles at many frequencies has 
been investigated. A clear trend has emerged which shows disagreement between the 
results from RPO and M-Layer at ranges near the horizon when surface based ducts are 
involved. In this thesis, results from three profiles are presented: a 300 m surface-based 
duct which is specified and used in Ref. 7 under the test name SBDUCT; a 14 meter 
evaporation duct which is also specified and used in Ref. 7, under the test name EDUCT; 
and a combination of these two ducts by merging the profile of the evaporation duct with 
that of the surface-based duct. The modified refractivity profile of the 300 m surface- 
based duct is given in Table m.l. The profile for the 14 m evaporation duct is given in 

Table n 12. The profile for the combination ^ ffl l A 300 m surface-based duct 
is given in Table IDL3. For each of the 
three profiles, three frequencies at 3 GHz, 

6 GHz and 12 GHz are selected. With the 
transmitter fixed at a height of IS m, the 
propagation loss of up to 100 m at ranges 
of 15, 20, 30, 40, 50, 60, 70, 80, 90, 100 
and 110 km is plotted. Since the results at different frequencies show similar features, 
only those at 6 GHz are included in this chapter. Those at 3 GHz and 12 GHz are 
collected in Appendices B through D. 


i 

Z| 

meters 


0 

0.000 

339.0 

1 

250.0 

368.5 

2 

300.0 

319.0 

3 

1000.0 

401.6 
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Table ni.2. A 14 m evaporation duct 


Zi 

meters 

Mi 

0.000 

339.00 

0.040 

335.10 

0.100 

333.66 

0.200 

332.60 

0.398 

331.54 

0.794 

330.51 

1.585 

329.53 

3.162 

328.65 

6.310 

327.96 

12.589 

327.68 

14.000 

327.67 

25.119 

328.13 

39.811 

329.25 

50.119 

330.18 

63.096 

331.44 

79.433 

333.12 

100.000 

335.33 

125.893 

338.20 

158.489 

341.92 

199.526 

346.69 

209.526 

347.87 


Table m.3. A 300 m surface-based 
duct over a 14 m evaporation duct 


Zi Mi 

meters 


0.000 339.00 


0.040 335.10 


0.100 333.66 


0.200 332.60 


0.398 331.54 


0.794 330.51 


1.585 329.53 


3.162 328.65 


6.310 327.96 


12.589 327.68 


14.000 327.67 


25.119 328.13 


39.811 329.25 


50.119 330.18 


63.096 331.44 


79.433 333.12 


100.000 335.33 


125.893 338.20 


158.489 341.92 


199.526 346.69 


209.526 347.87 


250.0 


300.0 319.0 


1000.0 401.6 


















































































































For all the eases, the polarization is chosen to be vertical. The receiver height 
increment delzr is set at 0.5 m. The receiver range increment delx is set at 2500 m. For 
RPO computations, the maximum height zmax is set at 100 m; the maximum range xmax 
is set at 115 km. The Gaussian beam pattern with a beamwidth of 5° is chosen for the 
antenna. Within the altitudes and ranges considered, there is no perceivable difference 
when the omni-directional pattern is used instead. For M-Layer computation, the 
parameter aloss is set to 2 dB/km. In all die figures, results from M-Layer are drawn as 
a solid line. Results from RPO are marked with asterisks if they are in the PE region and 
with dots if they are in the RO region. For easy reference, the height of the transmitter 
horizon at each distance, based on the four-thirds effective earth radius, is indicated with 
a horizontal line drawn across the figure on which it is present It starts to appear in the 
figure for the 20 km range. Figures beyond 50 km lie below the horizon completely and 
this horizontal line cannot be seen. 

In what follows, results of RPO and M-Layer computations at 6 GHz are presented. 
Their analysis and discussions are given in Chapter IV. 

A. 300 METER SURFACE-BASED DUCT 

Figures 2 through 12 show the propagation loss at various ranges from 15 through 
1 10 km when a 300 m surface-based duct is present This is the profile in which most 
significant deviations between RPO and M-Layer are observed, especially at 50 and 60 
km ranges. Even at 40 km, the two differ by 10 to 30 dB at low altitudes. On the other 
hand, the two agree well within line-of-sight and deep shadow regions. 
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Figure 4. Propagation loss at 30 km. 
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Figure 8. Propagation loss at 70 km. 



Figure 9. Propagation loss at 80 km. 



















Height (at) 



Figure 10. Propagation loss at 90 km. 



Figure 11. Propagation loss at 100 km. 
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Figure 12 . Propagation loss at 110 km. 


B. 14 METER EVAPORATION DUCT 

Figures 13 through 23 show the propagation loss in the presence of a 14 m 
evaporation duct RPO and M-Layer agree well over the entire range. From the figures 
in Appendix C , it can be seen that this is true for all the frequencies investigated, with 
only a less than 1.5 dB difference around where the loss is maximum for the 12 GHz case 
in the over-the-horizon region. 
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Figure 13. Propagation loss at 15 km. 
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Figure 15. Propagation loss at 30 km. 



Figure 16. Propagation loss at 40 km. 
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Figure 18. Propagation loss at 60 km. 
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21. Propagation loss at 90 km. 














Figure 23. Propagation loss at 110 km. 

C. SURFACE-BASED DUCT OVER EVAPORATION DUCT 

The agreement in one case and the disagreement in another between RPO and M- 
Layer prompted an investigation into a combined profile consisting of these two ducts. 
Figures 24 through 34 show the propagation loss under such a refiractivity profile. At 
6 GHz, RPO and M-Layer agree well over the entire range except between 50 and 60 km 
when M-Layer displays an increase in field strength toward the surface over the lowest 
few meters while RPO continues to decrease. This leads to more than 10 dB deviations 
in propagation loss at these ranges. At 12 GHz, the two programs agree even better when 
M-Layer does not show the increase toward surface level within the last few meters. At 
3 GHz, this difference is more pronounced. Below 10 to 20 m, this divergence starts to 
appear and leads to deviations in propagation loss of up to 40 dB. 
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Igure 25. Propagation loss at 20 km. 
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Figure 27. Propagation loss at 40 km. 
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Figure 29. Propagation loss at 60 km. 
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e 30. Propagation loss at 70 km. 



Figure 31. Propagation loss at 80 km. 
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IV. ANALYSIS AND CONCLUSIONS 


A. ANALYSIS 

Since RPO is a high frequency approximation while M-Layer is based on a low 
frequency technique, and the vertical extent within which the index of refraction of the 
atmosphere is specified is much larger than the wavelength under consideration, RPO is 
expected to be more accurate at short ranges. The fact that M-Layer matches the results 

m 

of ray-opdcs computations in Fig. 13 at 15 km from the transmitter for heights over 90 
m above ground confirms the reliability of M-Layer. RPO and M-Layer agree well also 
in the region far beyond the radar horizon. What happens then over the region in-between 
when a 300 m surface-based duct is present? 

The physics of the situation is revealed if the earth-flattening approximation as 
prescribed in Ref. 6 is re-examined. In fact, both RPO and M-Layer are treating die 
propagation problem as one involving only a flat "earth" surface. The waves are traveling 
in a flat, layered dielectric waveguide. There is no blockage of rays due to the earth to 
create a shadow; instead, the rays are bent away from the "shadow" region by the linearly 
increasing component of the modified refractivity. The ducting structure bends the wave 
back towards the surface. Properly implemented, these two programs should provide 
identical predictions over ranges of common validity. 

Above a flat surface, the concept of interference between a direct ray and a reflected 
ray, even though bent by the atmosphere, remains valid. The bending shifts the locations 
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of constructive and destructive interferences, but their separations are less affected. From 
Figs. 2 and. 13, it is clear that at 15 km, the separations between two neighboring nulls 
above ground are about 28 m, independent of the particular environment and close to a 
free-space, parallel-ray estimate of 25 m. The separations of nulls are expected to increase 
as the range is increased. This is observed in Figs. 13 through 23 for the propagation in 
the presence of the 14 m evaporation duct But when the 300 m surface-based duct is 
present, both RPO and M-Layer predict, for ranges greater than 90 km, similar field 
strength variations which oscillate much faster than those observed within the line-of- . 
sight region. It is clear that this is a phenomenon due to the presence of the surface-based 
duct In terms of high frequency ray-optics terminology, this fast oscillation is due to 
additional rays which are bent back by the duct In terms of waveguide mode theory, it 
is due to the propagation of waves in many modes established in the surface-based duct, 
each at a distinct phase velocity, thus interfering severely. 

The deviations between RPO and M-Layer at ranges from 30 through 70 km which 
show up in Figs. 4 to 12 can now be explained. RPO fails to fully take into account the 
effects of a high surface-based duct This argument is further supported by investigating 
the effects of lowering the duct height At all frequencies, the fast oscillation 
characteristic of the presence of a surface based duct sets in at a shorter range with a 
lower duct height Agreement between RPO and M-Layer gets better as the duct height 
is lowered. These can be observed from the results presented in Appendix E for 
propagation in a 150 m surface-based duct and in Appendix F for propagation in a 100 m 
duct 
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B. CONCLUSIONS 

At low altitudes in the neighborhood of the radar horizon, deviations between RPO 
and M-Layer can be significant. RPO appears not to be able to properly include the 
effects of a high surface-based duct at a short range. It is recommended that parabolic 
equation computations should start at a closer range than currently prescribed in RPO. 
The altitude at which the filtering of Held strength starts should always be much higher 
than the duct height 
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APPENDIX A: PROGRAM SOURCE CODES 
This Appendix contains listings of the RPO main program RPOmain, the RPO input 
subroutine RPOstdin, the M-Layer input subroutine MLstdin and the MATLAB M-file 
RPOMLA.M for plotting the propagation loss computed with RPO and M-Layer. 

1. Subroutine RPOmain 


program RPOmain 

c 

c PURPOSE: Main program for RPO. 

c 

c 

c input: 

c Argument List: None 

c Common: fmhz, ztran, ipolar, ipatrn, beamw, elang, hfang, hffac, nfacs 
c capm, zprof, nprofs, nlevls, xprof, wind 

c iscatt, maxlev, maxnx, nxout, maxnz, nzout, maxpro, xmax, zmax 

c 

c 

c output: 

c Argument list: None 
c Common: losscb, srng 
c 
c 

c subroutines calling RPOmain: None. 

c subroutines called by RPOmain: rpostdin, rpoini, rpostp, 

c 

c 

c common block 


c 

/system/ 

c 

fmhz : 

EM system frequency 

c 

ztran : 

antenna height 

c 

ipolar: 

antenna polarization 

c 

ipatrn: 

antenna pattern 

c 

beamw : 

antenna vertical beam width 

c 

elang : 

antenna elevation angle 

c 

hfang : 

height-finder angles array in degrees (0. to 99) 

c 

hffac : 

height-finder power reduction factor array (0. to 1.0 ) 

c 

nfacs : 

number of power reduction angles/factors for user-defined 

c 


height-finder radar 

c 

c 

/enviro/ 

c 

capm : 

(*) profile modified refractivity array 

c 

zprof : 

(*) profile heights in meters (.GE. 0.) 

c 

nprofs: 

number of profile levels 

c 

nlevls: 

number of refractivity profile levels (1 to maxlev) 

c 

xprof : 

(*) range to each profile in meter (.GE. 0.) 

c 

c 

wind : 

wind speed at range zero in meter/sec (.GE. 0.) 
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c 

/init/ 

c 

iscatt 

c 

maxlev 

c 

maxnx 

c 

nxout 

c 

maxnz 

c 

nzout 

c 

maxpro 

c 

xmax 

c 

zmax 


c Declares the 


flag to include troposcatter 

maximum number of profile level above zero height 
maximum number of output range points 
number of output range points (1 to maxnx) 
maximum number of output height points 
number of output height points (1 to maxnz) 
maximun number of profiles beyond zero range 
maximum range for output in meter 
maximum height for output in meter 
types of parameters 


character^ filein 


integer*4 ipolar, ipatrn, nfacs, nprofs, nlevls, iscatt, 
♦ nxout, nzout, maxlev, maxnx, maxnz, maxpro 

integer*2 losscb, srng 

real fmhz, ztran, beamw, elang, hfang, hffac, 

+ xprof, wind, capm, zprof, xmax, zmax 

c The following include file contains a PARAMETER statement to 
c define maximun array sizes maxpro, maxlev, maxnz 

INCLUDE 'RPOSIZE.INC' 


c 

c The following PARAMETER statement defines maximum array 
c dimensions used throughout RPO. Generally, these constants 

c will have to be changed for each implementation of RPO. 

c GLOSSARY: 

c maxlev: maximum number of profile levels above zero height, 

c maxnx : maximum number of output range points, 

c maxnz : maximum number of output height points, 

c maxpro: maximum number of profiles beyond zero range. 

c 

PARAMETER (maxnx = 440, maxnz = 280, maxlev = 50, maxpro =32) 
c 

dimension losscb(maxnz),srng(maxnz) 

COMMON /system/ fmhz, ztran, ipolar, ipatrn, beamw, elang, 

+ hfang(10), hffac(10), nfacs 

COMMON /enviro/ nlevls, nprofs, wind, zprof(0:maxlev, 0:maxpro), 
+ capro(0:maxlev, 0:maxpro), xprof(0:maxpro) 

COMMON /init/ iscatt, nxout, nzout, xmax, zmax 

COMMON /misc/ jminPE, nx, xstep, zstep, wl, rk, fterm, pi 

COMMON /inout/ filein, selx, mlxout 


c read data from input files 

write (*,*) 'Begin with rpostdin' 
call rpostdin 

c initialize RPO 

write (*,*) 'Begin with rpoini' 
call rpoini(nsteps) 


c call rpostp to compute propagation loss 
c and write the result to a file on disk 
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open(16,filesfilein//'.out') 
do iai, nxout 

call rpostp(x,losscb,srng) 

.Xkm*x*. 001 

write(*,*) 'Begin with rpostp',i.xkm 
if (x.ge.selx) than 
do jxl.nzout 

writa (16,1102) xkra, j*zstep, .l*losscb(j), srng(j) 
and do 
and if 
and do 

1102 format(f9.2, 6x. f6.1, 5x. f6.2. 5x, i2) 
and 


2. Subroutine RPOstdin 


c Subroutine RPOstdin raads in data* from a fila concocted from two 
c files: one contains parameters specific to RPO, the other contains 
c and parameters used by both RPO and M-Layer, including the modified 

c refractivity, the transmitter height, and the output points. 
c *************«****•**»**»*#***•************»****•*************•»#***! 

subroutine RPOstdin 

c The following include file contains a PARAMETER statement to 
c define maximun array sizes maxpro, maxlev 
INCLUDE 'RPOSXZE.INC' 
c 

c The following PARAMETER statement defines maximum array 
c dimensions used throughout RPO. Generally, these constants 

c will have to be changed for each implementation of RPO. 

c GLOSSARY: 

c maxlev: maximum number of profile levels above zero height, 

c maxnx : maximum number of output range points, 

c maxnz : maximum number of output height points, 

c maxpro: maximum number of profiles beyond zero range, 

c 

PARAMETER (maxnx = 440, maxnz = 280, maxlev = 50, maxpro =32) 
c 

character*8 filein, fileb 

COMMON /system/ fmhz, ztran, ipolar, ipatrn, beamw, elang, 

♦ hfang(10), hffac(10), nfacs 

COMMON /enviro/ nlevls, nprofs, wind, zprof(0:maxlev, 0:maxpro), 
+ capm(0:maxlev, 0:maxpro), xprof(0:maxpro) 

COMMON / in.it/ iscatt, nxout, nzout, xmax, zmax 
COMMON /inout/ filein, selx, mlxout 

c-read the RPO parameter- 

read(*,'(a)') filein 
read(*,*) nxout 
read(*,*) selx 
read(*, *) xprof(0) 
read(*,*) nprofs 

c RPO starts its array index with 0. It adopts Microsoft BASIC 
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c convention in its coding of this FORTRAN program. In the documentation 
c and the input file, "nprofs* is the total number of profiles. Thus 
c “nprofs* as the array index as used in RPO has to be adjusted by 
c substracting 'nprofs* by 1. 
nprofs*nprofs-l 
read(*,*) ipatrn 
read(*, *) beamw 

c For omnidirectional pattern, the parameter “beamw* is ignored in the 
c program. On the other hand, the upper limit for 'beamw* is 45 degrees, 
if (ipatrn.eq.l) then 
beamws45.0 
end if 

read( # ,*) elang 
read(*,*) nfacs 
read(*,*) iscatt 

c- read parameters common to RPO and M-Layer - 

read(*,*(a)') fileb 
read( # ,*) nlevls 

read(*, *) wind , 

read(*,*) fmhz 
read(*,*) ipolar 

c The definition of ipolar in RPO is 1 plus that in M-LAYER. ' 
ipolar=ipolar+l 
read(*,*) ztran 
read(*,*) nzout 
read(*,*) delzr 
read(*,*j zrinit 
zmax=nzout*delzr+zrinit 
read(*,*) mlxout 
read(*,*) delx 
read(*, *) xinit 
delx=delx*1000 
xinit=xinit*1000 
xmax=mlxout*delx+xinit 
do i^O,nlevls 

read(*,*) zprof(i,0) 
read(*,*) capm(i,0) 
capm(i,0)=capm(i,0)+339.0 
end do 

open(17,filesfilein//'.in') 

write(17,*)'-parameters common to RPO and M-Layer-' 

write(17,*)filein, ' :RPO file' 

write(17,*)nlevls, ' :nlevls (count from 0)' 

write(17,190)wind 
190 format(fl5.2, ' :wind') 
write(17,192)fmhz 
192 format(fl5.2, ' :fmhz') 

write(17,*)ipolar-1, ' :ipolar (0thorizontal l:vertical )' 

write(17,194)ztran 
194 format(f15.2, ' :ztran') 

writs(17,*)nzout, ' :nzout' 

write(17,201)delzr 

201 format(f15.2, ' :delzr') 
write(17,202)zrinit 

202 format(f15.2, ' :zrinit') 

write(17,*)mlxout ,' :mlxout (max number of range output) ' 

write(17,203)delx 

203 format(fl5.2, ' :delx') 





write(17,204)xinit 

204 forroat(£15.2, ' :xinit ‘) 

write (17, *) '- RPO parameters -' 

write(17,205)zmax 

205 format(£15.2, ' :zmax (nzout'delzr,calculated by program)' ) 
write(17,206)xmax 

206 format(fl5.2, ' :xroax (xinit+delx*mlxout,calculated by program)') 

write(17,*)nxout,' :nxout (xmax / delx, read from input) * 

write(17,207)selx 

207 format(f15.2, ' :selx (selected range x >= selx )') 

write(17,208)xprof(0) 

208 format(fl5.2, ' sxprof') 

write(17,*)nprofs+l, ' snprofs' 

write(17,*)ipatrn, ' sipatrn' 

write(17,209)beamw 

209 format(fl5.2, ' :beamw') 

write(17,210)elang 

210 format(f15.2, ' :elang') 

write(17,*) ' N/A :hfang (not used)’ 

write(17,*) ' N/A :hffac (not used)' 

write(17,*)nfacs, ' :nfacs' 

write(17,*)iscatt, ' :iscatt' 

write(17,*)'-M-profile-value on surface: 339.0 -' 

do i=0,nlevls 

write(17,211)zprof(i,0) 

211 format(f15.2, ' :zprof') 
write(17,212)capm(i,0) 

212 format(£15.2, ' :capm ') 
end do 

return 

end 


3. Subroutine MLstdin 


£*•*•*****•******•**••****•****•*•***•**'«••••*•*•**•*«•**•***•****•**•* 
c Subroutine MLstdin reads data from a file concocted from two files: * 
c one contains parameters specific to M-Layer, the other contains * 

c parameters used by both RPO and M-Layer, including the modified * 

c refractivity, the transmitter height, and the output points. * 

(.•IT******************************************************************** 

subroutine MLstdin 

c 

c MLstdin is the revised input program of M-Layer, NPS version, to 

c read in data files. The common block /inpt9/ has been removed from 

c all subroutines, 
c 

implicit real*8 (a-h,o-z) 
complex*16 qeigen 

integer iflgab.mpol,nzt,nzr,nx,nzlayr,nrmode.mfile,i 
character*40 filein, fileb 
c 

c use include file for parameters of 
c mxlayr max # layers 

c mxmode max # inodes 

c 

include 'mlaparm.inc' 
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c 

c 

c 

c 

c 

c 


c 


c 

c 

c 


c 


include file to define the 

maximum # of layers (mxlayr) 
maximum # of modes (mxmode) 


parameter (mxlayr=35 ) 
parameter (mxmode=390) 

dimension zi(mxlayr+1),zim(mxlayr*l),zigab(mxlayr*1), 
$ qeigen(mxmode) 


common 

♦ 

+• 

♦ 

* 

+ 


/inptO/filein.mfile,fqmzin.mpol,aloss,dielcg,sigmag, 

ztinit,delzt,nzt,zrinit,delzr,nzr,xinit,delx,nx, 
refz,refm,refgab,zim,zigab 
/inptl/nzlayr 
/inpt2/zi 
/inpt4/rmsbht 
/modes/nrmode,qeigen 


read M-LAYER specific parameters 
read(*,'(a)') filein 
read(*,*) mfile 
read(*,*) delfq 
read(*,*) nfreq 
read(*,*) aloss 
read(*,*) iflgab 
read(*,*) delzt 
read(*,*) nzt 
read(*,*) zref 
read(*,*) refz 
read(*,*) refm 
read(*,*) refgab 

read parameters common to RPO and M-Layer 
read(*,'(a)') fileb 
read(*,*) nzlayr 
read(*,*) wind 
read(*,*) fqmzin 
read(*,*) mpol 
read(*,*) ztinit 
read(*,*) nzr 
read(*,*) delzr 
read(*,*) zrinit 
read(*,*) nx 
read(*,*) delx 
read!*,*) xinit 


c The profile must contain at least three levels. The M gradients 

c in adjacent layers must not be equal, 

c 

do i*l,nzlayr+l 
read(*,*) zi(i) 
read{*,*) zim(i) 
end do 

if (mfile.ne.O) then 

c mfileoO indicates input file contains eigenvalues 
read(*,*)nrmode 
do i«l,nrmode 
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read(*,*) qeigen(i) 

•nd do 
endi£ 

c Calculate the root mean square bump height of sea surface. 

c 

rmsbhtsO.00514*wind**2 
c 

c Calculate the dielectric constant and conductivity of sea water. 

if (fqmzin .LB. 1500.) then 
dielcg= 80. 
sigma 3 4.3 

else if (fqmzin .LB. 3000.) then 

dielcgs 80. - .00733 * (fqmzin-1500.) 
sigma * 4.3 + .00148 * (fqmzin-1500.) 
else if (fqmzin .LB. 10000.) then 

dielcgs 69. - .00243 * (fqmzin-3000.) 
sigma * 6.52 + .001314 * (fqmzin-3000.) 
else 

dielcgs 51.99 
sigma 3 15.718 
end if 

c-write input data to a file on the disk- 

open(27,filesfilein//'.in') 

write(27,*) '- parameters common to RPO and M-Layer - 

write(27,' (2(A13))')filein. ':ML file' 
write(27,*)nzlayr, ' snzlayr' 

write(27,190)wind 
190 format(f15.2, ' :wind') 

write(27,192) fcpazin 
192 format(f15.2, ' :fqmzin') 

write(27,*)mpol, ' :mpol (0:horizontal lrvertical )' 

write{27,194)ztinit 
194 format(fl5.2, ' :ztinit') 

write(27,*)nzr, ' snzr' 

write(27,201)delzr 

201 format(£15.2, ' :delzr') 

write(27,202)zrinit 

202 format(fl5.2, ' :zrinit') 

write(27,*)nx, ' :nx' 

write(27,203)delx 

203 format(fl5.2, ' :delx in km') 

write(27,204)xinit 

204 format (f 15.2, ' .-xinit in km') 

write(27,*) '- M-Layer specific parameters - 

write(27,*) mfile, ' :mfilesOiread input and compute e.v.' 

write(27,205)delfq 

205 format(£15.2, ' :delfq') 

write(27,*)nfreq,' :nfreq' 

write(27,207)aloss 

207 format(fl5.2, ' :aloss') 

write(27,*)iflgab,' siflgab' 

write(27,208)delzt 

208 format(f15.2, ' :delzt') 

write{27,*)nzt, ' :nzt' 
write(27,209)zref 

209 format(fl5.2, ' :zref') 
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210 

211 

212 

214 

215 


write(27,210)refz 
forraat(fl5.2, ' :ra£z') 

write(27,211)refm 
format(£15.2, ' :ra£m') 

writa(27.212) refgab 
format(£15.2, ' :refgab') 

writa(27, *) '-M-pro£ile- 

do i*l,nzlayr+1 

write(27,214) zi(i),i 
format(£15.2, ' :zi(' ,12, ')' ) 

write(27,215) zim(i).i 
format(£15.2, ' :zim(' ,12, ')' ) 
end do 
return 
end 


4. MATLAB Plotting M-file 


% file name:"RPOMLA.M 
clg 

clear l 

fno=input(' input the no. of test file: ','s'); 

% delete the old .met file 

xdel*['delete rl', fno, '.met']; 
eval(xdel); 

% load the data 

tl=['load d:\matlab\wu\dat\', 'rpo', fno, '.dat'] 
t2=('load d:\matlab\wu\dat\','ml', fno, '.dat'); 
eval(tl); 
eval(t2); 

% read the data 

rpfr*['rpo', fno]; 
mlfr»['ml', fno]; 

rl=(rpfr,'(:,1)']; r2«(rpfr,'(:,2)']; 

r3=(rpfr,'{:,3)'] ; r4=(rpfr,'(:,4)'J; 
mlst(mlfr, ' (:, 1) ' ]; m2=(mlfr, ' (:, 3) ' J; m3s(mlfr, ' (:, 6) ' ] ; 
rpoxseval(r1); rpozaeval(r2); rpos=eval(r3); rporaeval(r4); 
mlx=eval(ml); mlz=eval(m2); mls=eval(m3); 

% set up initial conditoin for vdiile loop 
ansi; 

prof*input('The name for profile: ','s'); 
frq=input('The frequency (in MHz): '); 
anhsinput('The antenna height (m): '}; 

% loop for range increment 
while an>0 

xrg* input v ' The x range ' ); 

% calculate the receiver antenna height of horizontal distance 
ae*6371000; 

anh2* ( xrg*1000-sqrt(2*ae*anh) )~2 /(2*ae); 
kl»l; k2*l; )c3 = 1; k4sl; p*l; 

% xl, x2... is used to decide which region should be plotted in plot. 
Xls0; x2=0; x3=0; x4=0; 

% find the index of rpo loss for x range equal to xrg 
tmpxrg*xrg+.05; 

xindex*£ind( (rpox>=xrg) & (rpox<=tmpxrg) ); 
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stri=xindex(l); endi=xindex( length{xindex) ); 

% find the min. end max. of rpo loss for x range equal to xrg 
xrppmin=min(rpos(xindex)); 
xrpomaxsmax(rpos(xindex)) ; 

% get rpo loss, receiver height, according to the method of computation 
% (fe, ro, xo & pe) 

for i=stri:endi 

if rpor(i)==1 

rpzl(kl)= rpoz(i); rpsl(kl)= rpos(i); rprl(kl)= rpor(i); 
klakl+1; 

Xlsl; 

end 

if "-por(i) ==2 

rpz2(k2)= rpoz(i); rps2(k2)= rpos(i); 

k2=k2+l; 

x2=2; 

end 

if rpor(i)==3 

rpz3(k3)= rpoz(i); rps3(k3)= rpos(i); 
k3=k3+l; 

X3=3; . 

end 

if rpor(i)==4 

rpz4(k4)= rpoz(i); rps4(k4)= rpos(i); 

k4=k4+l; 

x4=4; 

end 

end 

% find the index of mlayer loss for the x range equal to xrg 
xindex=find( (mlx>=xrg) & (mlx<=tmpxrg) ); 
stri=xindex(l); endi=xindex( length(xindex) ); 

% find the min. and max. of ml loss for x range is equal to xrg 
xmlmin=min(mis(xindex)); 
xmlmax=max(mis(xindex)); 

% get the data which satisfy x=xrng 
for is3tri:endi 

mx(p)= mlx(i); mz(p)= mlz(i); ms(p)= mls(i); 

P=P+1; 

end 

% set the «ucis of plot 
raaxz=max(mz); 
maxx=max(xmlmax,xrpomax); 
minx=min(xmlmin,xrpomin); 
axis( [ minx maxx 0 maxz ] ) 
plot(ms,mz,'-g') 
hold on 

% plot the horizontal line 
hrs=(];hrz=(]; 
hrs=minx:1:maxx; 
for i*l:length(hrs) 
hrz(i)=anh2; 

end 

plot(hrs,hrz) 
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% plot rpo 

if xl*=l 

plot(rpsl,rpzl,'+b') 
end- 

if x2==2 

plot(rps2,rpz2,' .w') 

end 

if x3 ==3 

plot( rps3,rpz3,'xb') 

end 

if x4ss4 

plot( rps4,rpz4,'*r') 

end 

xlabel(' Propagation Loss (dB) '); 
ylabel(' Height (m) ') ; 
gtext([num2str(prof)]) ; 

gtext ([ '<-RPO' J) ; 

gtext ([' <-MLAYER' ]) ; 

gtext(['range : ' , num2str(xrg),'km']); 
gtext(['frequency: ',num2str(frq),'MHz']); 
gtext(['TX height: ',num2str(anh),'m']); 
if xl==l 

gtext(['FE region: +']) 
end 

if x2==2 

gtext(['RO region: . ‘ ) 
end 

if x3==3 

gtext(['XO region: x']) 
end 

if x4==4 

gtext(('PE region: *' ]) 
end 
grid 
pause 

csav=input(' save the plot ? (y/n) ','s'); 
if (csav=='y')|(csav=='Y') 
t3=('meta ', 'rl'.fno]; 
eval(t3); 

end 

ans=input(' enter N to Exit ','3'); 
if (ans=='N')i( ans=='n') 
an=0; 

end 

clear mz ms; 

clear rpsl rps2 rps3 rps4 
clear rpzl rpz2 rpz3 rpz4 
clg 

hold off 

end 
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APPENDIX B: PROPAGATION LOSS UNDER THE INFLUENCE OF A 300 M 

SURFACE-BASED DUCT 

This Appendix displays the propagation loss computed by RPO and M-Layer under 
the influence of a 300 m surface-based duct at 3 GHz and 12 GHz at ranges of IS, 20, 
30, 40, 50, 60, 70, 80, 90, 100 and 110 km. 

1. Propagation loss at 3 GHz 

Figures B.l through B.l 1 displays the propagation loss at 3 GHz computed by 
RPO and M-Layer. 
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Figure BJ. Propagation loss at 20 km. 
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Figure BJ. Propagation loss at 30 km. 



Figure B.4. Propagation loss at 40 km. 
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Figure B.5. Propagation loss at SO km. 
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Figure B.6. Propagation loss at 60 km. 
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Figure B.8. Propagation loss at 80 km. 
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Figure B.ll. Propagation loss at 110 km. 


2. Propagation loss at 12 GHz 

Figures B.12 through B.22 displays the propagation loss at 12 GHz computed 
by RPO and M-Layer. 


54 



























Height (m) 


300m Surface-Based Duct 



Figure B.14. Propagation loss at 30 km. 
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Figure B 20 . Propagation loss at 90 km. 



Figure B.21. Propagation loss at 100 km. 
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Figure B 22. Propagation loss at 110 km. 
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APPENDIX C: PROPAGATION LOSS UNDER THE INFLUENCE OF A 14 M 

EVAPORATION DUCT 

This Appendix displays the propagation loss computed by RPO and M-Layer under 
the influence of a 14 m evaporation duct at 3 GHz and 12 GHz at ranges of 15, 20, 30, 
40, 50, 60, 70, 80, 90, 100 and 110 km. 

1. Propagation loss at 3 GHz 

Figures C.l through C.11 displays the propagation loss at 3 GHz computed by 
RPO and M-Layer. 
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Figure C.l. Propagation loss at 15 km. 
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Figure C 2, Propagation loss at 20 km. 
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ire C J. Propagation loss at 30 km. 
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Figure C.4. Propagation loss at 40 km. 
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Figure C j 6 . Propagation loss at 60 km. 
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2. Propagation loss at 12 GHz 

Figures C.12 through C.22 displays the propagation loss at 12 GHz computed 
by RPO and M-Laycr. 
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Figure C.13. Propagation loss at 20 km. 
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re C.16. Propagation loss at SO km. 



Figure C.17. Propagation loss at 60 km. 
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C.18. Propagation loss at 70 km. 
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Figure C.19. Propagation loss at 80 km. 
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APPENDIX D: PROPAGATION LOSS UNDER THE INFLUENCE OF A 300 M 


SURFACE-BASED DUCT OVER A 14 M EVAPORATION DUCT 
This Appendix displays the propagation loss computed by RPO and M-Layer under 
the influence of a 300 m surface-based duct over a 14 m evaporation duct at 3 GHz and 
12 GHz at ranges of 15, 20, 30, 40, 50, 60, 70, 80, 90. 100 and 110 km. 

1. Propagation loss at 3 GHz 

m 

Figures D.l through D.l 1 displays the propagation loss at 3 GHz computed by 
RPO and M-Layer. 
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Figure D.2. Propagation loss at 20 km. 
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Figure DJ. Propagation loss at 30 km. 
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Figure D.4. Propagation loss at 40 km. 
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Figure D.6. Propagation loss at 60 km. 
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Figure DJJ. Propagation loss at 80 km. 
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Figure D.9. Propagation loss at 90 km. 



Figure D.10. Propagation loss at 100 km. 
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Figure D.ll. Propagation loss at 110 km. 


2. Propagation loss at 12 GHz 

Figures D.12 through D.22 displays the propagation loss at 12 GHz computed 
by RPO and M-Layer. 
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Figure D.15. Propagation loss at 40 km. 
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Figure D.16. Propagation loss at 50 km. 
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Figure D.18. Propagation loss at 70 km. 



Figure D.19. Propagation loss at 80 km. 
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Figure D 22 . Propagation loss at 110 km. 
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APPENDIX E: PROPAGATION LOSS UNDER THE INFLUENCE OF A 150 M 

SURFACE-BASED DUCT 

This Appendix displays the propagation loss computed by RPO and M-Layer under 
the influence of a 150 m surface-based duct at 3 GHz, 6 GHz and 12 GHz at ranges of 
15, 20, 25, 30, 35, 40, 50, 60, 70, 80, 90, 100 and 110 km. The modified refractivity 
profile is given in Table E.l. 


Table E.l. A 150 m surface-based duct. 


1 

Z 

meters 

Mi 

, 

0 

0.0 

339.0 

1 

100.0 

350.8 

2 

150.0 

301.3 

3 

1,000 

401.6 


1. Propagation loss at 3 GHz 

Figures E.1 through E.13 displays the propagation loss at 3 GHz computed by 
RPO and M-Layer. 
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Figure E.5. Propagation loss at 35 km. 
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Figure EUJ. Propagation loss at 40 km. 
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E.11. Propagation loss at 90 km. 
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Figure E.12. Propagation loss at 100 km. 
























Figure E.13. Propagation loss at 110 km. 


2. Propagation loss at 6 GHz 

Figures E. 14 through E.26 displays the propagation loss at 6 GHz computed 
by RPO and M-Layer. 


94 































Height (iu) 



Figure E.16. Propagation loss at 25 km. 
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Figure E.17. Propagation loss at 30 km. 
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re E.18. Propagation loss at 35 km. 



Figure E.19. Propagation loss at 40 km. 


97 































































Propagation Lon (dB) 
Figure EL23. Propagation loss at 80 Ion. 
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3. Propagation loss at 12 GHz 

Figures E.27 through E.39 displays the propagation loss at 12 GHz computed 
by RPO and M-Layer. 


101 


























Height (in 



15(hn Surface-Based Duct 



135 140 145 150 155 160 165 170 

Propagation Loss (dB) 

Figure 1L28. Propagation loss at 20 km. 
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Figure EJ4. Propagation loss at 60 km. 




























































Figure EJ5. Propagation loss at 70 km. 
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Figure EL36. Propagation loss at 80 km. 
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Figure EJ8. Propagation loss at 100 km. 
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Figure EJ9. Propagation loss at 110 km. 
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APPENDIX F: PROPAGATION LOSS UNDER THE INFLUENCE OF A 100 M 

SURFACE-BASED DUCT 

This Appendix displays the propagation loss computed by RPO and M-Layer under 
the influence of a 100 m surface-based duct at 3 GHz, 6 GHz and 12 GHz at ranges of 
15, 20, 25, 30, 35, 40, 50, 60, 70, 80, 90, 100 and 110 km. The modified refractivity 
profile is given in Table F.l. 


Table F.l. A 100 m surface-based duct 


i 

Zi 

meters 

M i I 

0 

0.0 

339.0 | 

1 

50.0 

344.9 

1 2 

100.0 

295.4 

1 3 

1000.0 

401.6 


1. Propagation loss at 3 GHz 

Figures F.l through F.13 displays the propagation loss at 3 GHz computed by 
RPO and M-Layer. 
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Figure FJ. Propagation loss at 35 km. 



Figure F.6. Propagation loss at 40 km. 
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Figure F.8. Propagation loss at 60 km. 
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ire F.9. Propagation loss at 70 km. 



Figure F.10. Propagation loss at 80 km. 
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Figure F.12. Propagation loss at 100 km. 




















2. Propagation loss at 6 GHz 

Figures F.14 through F.26 displays the propagation loss at 6 GHz computed 
by RPO and M-Layer. 
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Figure F.19. Propagation loss at 40 km. 
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Figure F.23. Propagation loss at 80 km. 
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Figure F26. Propagation loss at 110 km. 


3. Propagation loss at 12 GHz 

Figures F.27 through F.39 displays the propagation loss at 12 GHz computed 
by RPO and M-Layer. 
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Figure F.31. Propagation loss at 35 km. 



126 









































Height (m) 
































Propagation Loss (dB) 


Figure F.35. Propagation loss at 70 km. 



Propagation Loss (dB) 


Figure F.36. Propagation loss at 80 km. 









































Figure F 38. Propagation loss at 100 km. 
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Figure F.39. Propagation loss at 110 km. 

























LIST OF REFERENCES 


1. H.V. Hitney and J.H. Richter, "Integrated refraction effects prediction system 
(IREPS),” Naval Engineers Journal, 257-262, 1976. 

2. W.L. Patterson, et al., "Engineer’s Refractive Effects Prediction System (EREPS) 
Revision 2.0," Technical Document 1342 , revision 2.0, Naval Ocean Systems 
Center, San Diego, CA 92152-5000, February 1990. 

3. H.V. Hitney, "Engineer’s refractive effects prediction system (EREPS), Discussion 
Summary," in "Proceedings: conference on microwave propagation in the marine 
boundary layer," Technical Report TR 89-02, p. 3-71, Naval Environmental 
Prediction Research Facility, Monterey, CA 93943-5006, January 1989. 

4. L.W. Yeoh, "An analysis of M-Layer a multilayer tropospheric propagation 
program," Technical Report NPS-62-90-009, Naval Postgraduate School, Monterey, 
CA 93943, June 1990. 

5. H.-M. Lee and Y.Y. Han, "M-Layer NPS Version," IEEE Transactions on 
Magnetics, 29(2), 1363-1367, 1993. 

6. D.E. Kerr, ed.. Propagation of Short Radio Waves, Peregrinus, London, 1987; 
original edition 1951. 

7. W.L. Patterson and H.V. Hitney, "Radio physical optics CSCI software documents,” 
Technical Document 2403, Naval Command, Control and Ocean Surveillance 
Center, RDT&E Division, San Diego, CA 92152, December 1992. 

8. F.D. Tappert, "The parabolic approximation method," in J.B. Keller and J.S. 
Papadakis, ed.. Wave Propagation and Underwater Acoustics, 224-287, Springer- 
Verlag, 1977. 

9. J.R. Kuttler and G.D. Dockery, "Theoretical description of the parabolic 
approximation/Fourier split-step method of representing electromagnetic propagation 
in the troposphere," Radio Science, 26(2), 381-393, 1991. 


131 





INITIAL DISTRIBUTION LIST 


1. Defense Technical Information Center 2 

Cameron Station 

Alexandria, VA 22304-6145 

2. Dudley Knox Library, Code 52 2 

Naval Postgraduate School 

Monterey, CA 93943-5002 

3. Chairman, Code EC 1 

Department of Electrical and Computer Engineering 

Naval Postgraduate School 
Monterey. CA 93943-5121 

4. Professor Kenneth L. Davidson (MR/Ds) 1 

Department of Meteorology 

Naval Postgraduate School 
Monterey, CA 93943-5100 

5. Professor Lawrence J. Ziomek (EC/Zm) 1 

Department of Electrical and Computer Engineering 

Naval Postgraduate School 
Monterey, CA 93943-5121 

6. Professor Hung-Mou Lee (EC/Lh) 2 

Department of Electrical and Computer Engineering 

Naval Postgraduate School 
Monterey, CA 93943-5121 

7. Mr. Yeoh, Lean Weng 1 

833 Dyer Road, Room 537C 

Department of Electrical and Computer Engineering 
Naval Postgraduate School 
Monterey, CA 93943-5121 

8. Mr. Wu, Chi-Wei 1 

No. 22, Lane 35, Wu-Hwa Street, San-Chung City, 

Taipei County, Taiwan, R.O.C 


132 



1 


9. Mr. Ting, Chueh 
Department of Electrical Engineering 
Chung Cheng Institute of Technology 
Tao-Yuan, Tai-Hsi 
Taiwan, R.O.C. 

10. Program Executive Office, Theater Air Defense (D-21) 1 

Atm: Mr. George Hamilton, NC 2, Room UN06 

2531 Jefferson Davis Highway 
Arlington, VA 22242-5170 

11. NATO Seasparrow Surface Missile System Project Office (N-US) 1 

Atm: LCDR Scott J. Smith, 4 Crytal Park, Room 210 

2531 Jefferson Davis Highway 
.Arlington, VA 22242-5170 

12. Naval Command, Control and Ocean Surveillance Center 1 

RDT&E Division 

Ocean and Atmospheric Sciences (Code 54) 

San Diego, CA 92152-5000 


133 


